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Abstract 

' The bi-Hamiltonian structure is established for the perturbation 

• equations of KdV hierarchy and thus the perturbation equations them- 

O ! selves provide also examples among typical soliton equations. Besides, 

I a more general bi-Hamiltonian integrable hierarchy is proposed and 

\^ ■ a remark is given for a generalization of the resulting perturbation 

0^ , equations to 1 + 2 dimensions. 

Perturbation theory is the study of the effects of small disturbances. Its 
basic idea is to find approximate solutions to a concrete problem by exploiting 
O ■ the presence of a small dimesionless parameter. There have been a lot of 

works to investigate the perturbated soliton equations (see for example 
PI and references therein). Tamizhmani and Lakshmanan have considered a 
^ ! perturbation effect of the unperturbated KdV equation and they have given 

I rise to infinitely many Lie-Backlund symmetries and a Hamiltonian structure 

for the resulting equations in Ref. However they haven't obtained the 
corresponding bi-Hamiltonian formulation. 

In this letter, we would like to discuss the perturbation equations of the 
unperturbated whole KdV hierarchy, i.e. the effect of the disturbance around 
solutions of the following original KdV hierarchy 

ut„ = Hu) =d^ + 2u,d-' + Au, = -^, n>l. (1) 

The first equation is exactly the usual KdV equation 

Uh = Uxxx + QuUx, (2) 
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which describes the unidirectional propagation of long waves of small ampli- 
tude and has a broad of applications in a number of physical contexts such as 
hydromagnetic waves, stratified internal waves, ion-acoustic waves, plasma 
physics and lattice dynamics (for details see, for example, 0). The op- 
erator is a common hereditary recursion operator for the whole KdV 
hierarchy. We shall show that all of the resulting perturbation equations pos- 
sess a bi-Hamiltonian structure and thus they constitute a typical integrable 
soliton hierarchy. We shall also point out a mistake on the Hamiltonian for- 
mulation in Ref . . Finally, a more general bi-Hamiltonian integrable soliton 
hierarchy is established and some further discussion is presented. 

We first consider the case of KdV equation and then consider the 
case of higher order KdV equations. Let us make a perturbation expansion 



N 



U 



1=0 



iV > 1, 



and call the equation 



oie 



N\ 



+ QuUx + o(e 



(3) 



(4) 



the A^-th order perturbation equation of KdV equation (|^). It is easy to find 
that 

N 
j=0 

where 5,0 is the Kronecker's symbol. Therefore the iV-th order perturbation 
equation (H) may be rewritten as 



Ut 



N N 



i=0 



i=0 



A balance of the coefficients of like powers of e leads to the following equiv- 
alent equation 

% = ^i.v)'nx, V = ivo, ■ ■ ■ , Vn)'^, (5) 

where the operator reads as 



$ = $(77) 



$(?7 



N 







$(r/i) ^(r/o) 



(6) 



By its first component of (^), we see that the perturbation expansion (^) is 
an expansion around an exact solution r/o of KdV equation (||). Moreover u 
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is an approximate solution of KdV equation (0) to a precision o(£:^) when rj 
satisfies (|^). However our aim here is not to get some approximate solutions. 
What we want to discuss is some algebraic structures that the perturbation 
equations possess. 

Let us introduce two differential operators 
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^(^i) 



N) 
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, J = 










_ 



(7) 



with J{r]i) = 6iod^ + 2rjix + irjid, <i < N. Note that the operator M + aJ 
is a Hamiltonian operator for any constant a (see for example Ref. [§) and 
that 

l>(r/) = MJ-\ = <lj, ^ = $*, 

where the asterisk appended to $ denotes the conjugate operation. There- 
fore J,M constitute a Hamiltonian pair and further the operator is 
hereditary 

We are now in the position to construct a bi-Hamiltonian structure |^ 



T0| for the perturbation equation (Rf). We first have the first Hamiltonian 



structure 



5Ho 
5rj 



M- 



5rj 



N 



Ho = Hoiv) = 7;J2ViVN^i- 



1=0 



Second noting that 

^(ri,) := ^*{r]i) = 5iod^ + 2T]i + 2d-%d, 0<i<N. 
we can compute 



/i(r/) := ^{ri) 



5rj 



N 







VNxx + 3 X] i + j = N f]if]j 

i,j>0 

W-l,a;a; + 3X] i + j = N-\ TJirjj 

i,j > 



VOxx + ^Vo 



6Hi 





Vn 




Vn-1 








^0 



6ri 
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where the Hamiltonian function Hi is defined by 
/■I 1 

*^ i + j = N i + j + k = N 

i,j>0 i,j,k>0 

Here and hereafter < -, ■ > stands for the standard inner on R'^^^. The above 
analysis allows us to conclude that the perturbation equation (j^) possesses 
a bi-Hamiltonian structure 



5Hi 
5rj 



M 



5rj 



(9) 



When iV = 1, N = 2, the perturbation equation (R) becomes 



^0 

ijietaoVox 
Vixxx + Q{'noVi)x 







Voxxx + QvoVox 






Vlxxx + G{VoVi)x 


. V2 _ 


ti 


_ Tl2xxx + ^VlVlx + G{V0V2)x _ 



respectively. These two equations are all integrable coupling ones with KdV 
equation. Integrable coupling is an interesting subject in soliton theory |1T 



Note that the second equation may be reduced to the first one if we take 
Vi = 0- 

For the case of higher order KdV equations, by induction we can engender 
the following equivalent equation 

Vt^ = ^{v)Vtr.-i = ^"{v)Vx (n > 1) (10) 
to the A^-th order perturbation equation of the n-th KdV equation 



7i,„ = $"(^)n.. + o(£^) 



'11^ 



namely, 



N N 
^[J2'^iV^y){T.V^t.^^^') (mode 



i=0 



i=0 



According to the bi-Hamiltonian theory (for a detailed description see |T2 
l3| and [Q), we know the vectors ^"(r^)^^, n > 1, are all gradient vectors 



because ^(r/)^^ is a gradient vector and ^* = $ is hereditary. Therefore 
there exist functions Hn so that 
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Moreover the Hamiltonian functions can be computed by 



Hn = Hn{ri) = f < fn{Xv),V> dX. 
Jo 



Therefore we see that the perturbation equation (|10D of the n-th KdV equa- 
tion in the hierarchy (|1|) possesses a bi-Hamiltonian structure 



6H„ 



M 



6K 



n-l 



'121 
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In fact, the perturbation equation (|10D has a multi-Hamiltonian structure 



Vt„ = J- 



5rj 



6ri 



. 5 Ho 
61] 



where J\l/\ < i < n, are all Hamiltonian operators and constitute Hamil- 
tonian pairs with each other. It is also interesting to note the hereditary 
operator $"(17) has the following concrete form 



A, 



A 



N 







Ai Ao 



where the operator Aj, < j < N, are determined by 

Aj= E ^V^^)■■■HV^J, 0<J <N. 

«i H \-in= j 



Unlike the hereditary operator for the Kepler system [|T5[, the above heredi- 
tary operator can generate new independent constants of motion Hn from a 
starting one Ho- 

By now, we have proposed a bi-Hamiltonian structure for the perturba- 
tion equations of the whole KdV hierarchy. This bi-Hamiltonian structure 
yields infinitely many symmetries and infinitely many corresponding con- 
stants of motion in involution for every perturbation equation and thus the 
perturbation equations of KdV hierarchy provide a typical hierarchy of soli- 
ton equations. It is worth noting that the algebraic structures of the first 
order perturbation equations may also be well described by the perturbation 
bundle and cotangent bundle fl^ . 



We remark that some cases, especially = 1, for KdV equation (0) has 
in detail been considered in Ref.|Q|. But they only described a Hamiltonian 
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structure. Moreover that Hamiltonian structure is not suitable for the per- 
turbation equations of KdV equation, which leads to some confusion of the 
symbols. For example, in the case of = 1, we cannot write 



d 
d 



Actually, there doesn't exist this kind of Hamiltonian function Hq, because 
iVoxx + ^v'^yVixx + GvoVi)'^ is not a gradient vector. However this is a small 
problem, which is easily changed to the correct version mentioned above. 

There also exists a more general result than the bi-Hamiltonian structure 



(|T^). Let as, 1 < s < p, 6^, c^, 1 < s < q, dg, 1 < s < r, are arbitrary real 
constants, and kg, I < s < p, Ig I < s < q, m^, 1 < s < r, are distinct 
non-negative integers not greater than A^. We further choose that J = J and 



M 







Avo) 
Jivi) 



(13) 



^(^o) Avi) ■■■ Avn) 

where the operator J{rii), < i < N, are defined by 

^ p q r 

s=l s=l s=l 

+2rji, + 4r]id, 0<z<N. (14) 
Note that ( p!^ means we choose 

AVks) = O'sd'^ + 2r]ks,x + 4?7fc^<9, 1 < s < p, 
A'nis) = bsd^ + Cgd + 2r]i^^^ + Ar]i^d, 1 < s < g, 

JiVms) = dsd + 27]m„x + '^Vm.d, 1 < S < T, 

J(?7i) = 2riinc + Ariid, {/c^, l< s <p] h, I < s < q; m^, 1 < s < m}. 

The differential operators J and M are still a Hamiltonian pair and hence 
the operator ^{rj) = MJ~^ is a hereditary symmetry operator. This leads 
to the following hierarchy of integrable equations which possesses a more 
general bi-Hamiltonian structure 



Vtn = = J- 



5r] 5r] 
with the Hamiltonian functions defined by 



M — , n> 1 



(15) 



< {^''Vx){Xv),V> dX, n>0, 
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in which Hq and Hi read as 

1 ^ 

Ho = Ho = -Y,ViVN-i, 
1 ^ 

^1 = 2 51 ViVjxx 

s=l i^j = N-ks 

i,j > 

1 

s=l i + j = N-ls i + j = N -Is 

1 

s=l i ^ j = N - i + j + k = N 

i,j>0 i,j,k>0 

Evidently, (|T2|) is a special case of the hierarchy (|15|), one reduction with 
ki = and ai = 1, = 0, 1, bg = Cg = ds = 0, Vs. 

We point out that our method can also be applied to discussing the forced 



integrable systems [0. It is sufficient to further make the perturbation of 
value boundary conditions and then we can get new forced integrable systems. 
The resulting forced systems can be solved by the same approach as the one 
in the original forced integrable systems Jl^. However in the present letter, 
we haven't considered this kind of interesting problems due to the limited 
space. 

Finally, we suggest that it may be more fruitful to think of the pertur- 
bation expansion u as depending on two independent space variables, x and 
y = ex, i.e. 

N 

^ = J2 Vii^^ ^2, ■ ■ y = ^x. (16) 

This kind of expansions comes from the perturbation method of multiple 
scales and the variable y is called a slow variable. The perturbation 



method of multiple scales for time variable has been applied to the study of 
perturbated soliton equations (see for instance p|) and the study of solutions 
to soliton equations (see for instance JT^). Under the expansion (0), we have 



d _ d d 
dx dx dy 

Therefore from the iV-th order perturbation equation of KdV equation 
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we can similarly obtain the following equation in 1 + 2 dimensions 

' ^70*1 = VOxxx + QVoVox, 

Vlh = Vlxxx + ^VOxxy + Q{V0lll)x + QVoVOy, 
< V2h = V2XXX + ^Ixxy + ^VOxyy + Q{VoV2)x + ^VlVlx + Q{VoVl)y, 
Vjti — Vjxxx + ^Vj-ljXxy + ^Vj-'ijXyy + Vj-S,yyy 

+6 (eLo ViVj-i,x + EiZo ViVj-i-i,y) , ^<j<N. 

This is a generalization of the perturbation equations of KdV hierarchy. How- 
ever, we don't know the complete answer on the integrability of the above 
equation, even in the case N — 1 where the above equation reads as 

VOh = VOxxx + QvoVox, 

Vlti = Vlxxx + ^Oxxy + Q{V0Vl)x + QVoVOy 

This equation is also a coupling one with KdV equation. If integrable, it 
provides us an another example which gives an integrable coupling equation 
with KdV equation. 
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